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(57) Abstract 

It is described a discrimination 
procedure of a wanted signal from a 
plurality of cochannel interferents re- 
ceived by array antennas of GSM or 
DCS base transceiver stations. The 
procedure includes a phase for the es- 
timate of the number and arrival di- 
rections of the interferents, and of the 
wanted signal, followed by a spatial fil- 
tering phase in which the signals trans- 
duced by the sensors of the array are 
linearly combined among them through 
multiplication coefficients, or weights, 
organized in a vector w satisfying the 
two following conditions: A) Spatial 
filtering constrains the gain of the array 
in the ratio between wanted signal and 
noise, compared to the traditional use 
of a single sensor, so that the gain is not 
less than a properly selected threshold; 
B) it minimizes the ratio between the 
sum of interferents' powers and wanted 
signal power. 
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DISCRIMINATION PROCEDURE OF A WANTED SIGNAL FROM A PLURALITY OF 
COCHANNEL INTERFERING SIGNALS AND RECEIVER USING THIS PROCEDURE 



Field of the invention 

The present invention relates to the field of the so-called "intelligent" array antennas, 
and more in particular to a discrimination procedure of a wanted signal from a plurality 
of cochannel interferents received by array antennas of base transceiver stations for 
0 cellular telecommunication and relative receiver. 
Background art 

The use in mobile radio environment of antenna consisting of one or more 
arrays of electromagnetic field sensors (array) is known, for instance, from the 
European patent application published under No. 0593822A1. This document in fact, 
describes, mentioning the first claim, "An arrangement of antennas for a base 
transceiver station including a plurality of antenna arrays, each array able to form a 
multiplicity of narrow radiation lobes, separated and partially overlapped in the 
azimuth plane, the arrays being positioned in such a way that all the iobes generated 
by the same give a substantia! omnidirectional coverage on said plane; lobes 
formation means (beamforming) in azimuth and angular elevation for each said array; 
a plurality of radiofrequency transceivers each one to transmit and receive the signals 
relative to one or more channels; a switching matrix to connect each transceiver with 
one or the other array through said beamforming means; etc.". 

More in particular, a beam former, for example a Butler matrix, is associated to 
an antenna array. A beam former consists of phase shifting and adder circuits, used 
in reciprocal way both in reception and in transmission. Considering for instance 
reception, the beam former has a plurality of input ports for the signals coming from 
the sensors of the array , and a multiplicity of output ports, each one relative to a pre- 
set azimuth direction and corresponding to a particular combination of module and 
phase of input signals. A dual behaviour applies to transmission, where the same port 
selected during reception is used to transmit towards the mobile . The antenna array, 
and relative beam former, are therefore an essential part of a system capable to 
identify the direction of the signal transmitted uplink by a mobile, that the system 
follows condensing a narrow radiation lobe in which the power of the signal 
transmitted down-link towards the mobile itself is concentrated. This equals, to an 



WO 99/33141 H H PCT/EP98/08165 



10 



J 5 



intelligent behaviour of the antenna, which deviates from the traditional utilization of 
antennas in the same sector of the technique. Thanks to tne intelligent behaviour, the 
interference from cochannel channel is reduced and the re-utilization of the same 
frequencies in adjacent cells is made possible. The deriving advantage is 
considerable and consists in the possibility to increase the dimension of cells, at equal 
transmitted power, in low traffic areas or reduce the frequency re-utilization distance, 
increasing the number of carriers per cell in high traffic areas. 

The estimate of wanted signal arrival direction avails of means for the 
measurement of signals present at output ports of the beam formers and of a 
processor evaluating the above mentioned measures and selecting the best direction, 
for instance that one for which the signal level is higher. 

The antenna arrangement described in the mentioned application, overcomes 
the limits introduced by sectorial antennas in corner-excited cells, still widely employed 
in the mobile radio field, due to need of infrasectorial handovers as the position of the 
mobile part around common antennas varies, though remaining close to the same. 
The intelligent antenna system enables in fact a transceiver to be associated to any 
possible narrow radiation lobes, sunburst arranged on the azimuth plane, therefore 
the cell appears to the network as if it were equipped of omnidirectional antenna, but 
without the relative known drawbacks. 
20 Background art 

Notwithstanding all the features mentioned above, the intelligent antenna 
arrangement described in the mentioned application, fails to exploit the capability of 
an antenna array to shape the overall radiation diagram so that nulls are steered in 
the interferents' directions. For. example, if interfering and wanted signal arrival 
directions are different but contained in the same beam, no discrimination can be 
made between them. This stands also if beam is selected, after demodulation, on the 
basis of wanted signal quality. 

In order to exploit the interference nulling capabilities of the antenna array a 
numerical (and not analog) beamforming has to be performed on digitalized signals 
30 from the sensors of the array. 

Numerical beamforming is nothing more then a linear combination of base band 
signals from the array elements using a set of complex coefficients w. 
The procedures for numerical beamforming mainly differ in the choice of the 
coefficient's set w and usually 2 different methods can be used (see [i] for a more 
detailed description and for a comparison based on system capacity): 
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1) All cochannel user arrival directions are estimated and coefficients w are chosen 
with the constraint that the overs!! array gain is 'nuiiect' in the interferent arrival 
directions. See also [2| for an overview of known art about beamforming and "null 
steering" 

5 2) Coefficients w are calculated so that the difference between the combined signal 
and a reference part of the signal is minimized. 
Both these methods for coefficients w calculation suffer from the fact that they 
optimize uplink receiver performance but the coefficients w are not always reliable for 
downlink transmission. 

10 For example, when uplink noise level is negligible and the interferent and wanted 
signal arrival directions are close each other, both methods try to reduce interference 
as much as possible. The result is often a very low level of the array gain in the 
direction of the wanted signal and thus only a small fraction of transmitted energy will 
be directed towards the user. 

15 Summary of the invention 

The scope of the present invention is to overcome the above mentioned 
drawbacks and to indicate a discrimination procedure of a wanted signal from a 
plurality of cochannel interferents received by array antennas of base transceiver 
stations for cellular telecommunication. 

20 To attain these objects, scope of the present invention is a discrimination 

procedure of a wanted signal from a plurality of cochannel interferents received by 
array antennas of frequency division (FDMA), or time division (TDMA), or mixed 
FDMA/TDMA multiple access telecommunication systems, re-employing a same 
frequency group in adjacent territorial areas, including an estimate phase of the arrival 

25 directions of the said interferents and of said wanted signal, and a successive phase 
of spatial filtering in which signals transduced by the relative sensors of a said array 
are linearly combined among them through multiplication coefficients, or weights, 
obtaining a reception signal cleaned from interferents; said phases being repeated for 
each one of the time slots, in the TDMA frame, in which the wanted user transmits, 

30 characterized in that said weights satisfy the two following conditions: 

A) the gain of said array in the ratio between said wanted signal and the noise after 
spatial filtering, compared to the traditional use of a single sensor, is constrained so 
that it is is not less than a properly selected threshold; 

B) the ratio between the sum of interferents* powers and wanted signal power is 
35 minimized, as described in claim 1. 
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The subject procedure can find useful application in the realization of a receiver 
for base transceiver stations of cellular telephone systems of the GSM 900 MHz, or 
DCS 1S00 MHz type. 

The occurrence of condition A) enables to obtain some advantages: 
5 1. to improve the performance in uplink reducing the performance sensitivity of the 
spatial filtering to the error made in the estimate of wanted signal arrival directions. 
2, To enable, in downlink, a high directivity of the antenna array radiation diagram in 
the direction of the wanted. 

The occurrence of condition B) enables to improve the performance in the two 
10 connections, uplink and downlink, reducing the interference. 

Profitably, it is possible to express in the vectorial form the equations 
concerning the estimate phases of the arrival direction and spatial filtering, as well as 
those generated by the imposition of the above mentioned featuring conditions, thus 
making the mathematical computation for the determination of said weights more 
15 compact. 

Further object of the present invention is a frequency division multiple access 
receiver for telecommunication systems (FDMA), or time division multiple access 
(TDMA), or mixed (FDMA/TDMA), which re-employ a same frequency group in 
adjacent territorial areas, characterized in that it includes means for the actuation of 

20 the discrimination procedure of a wanted signal from a plurality of cochannel 
interferents, which already formed the object of an invention, as described in claim 13. 
Brief description of drawings 

Additional scopes and advantages of the present invention will be more evident 
from the following detailed description of an embodiment of the same and attached 

25 drawings, in which: 

- Figures 1 and 2 show an array of sensors used in the receiver of the present 
invention, when it is invested by plane waves coming from different directions; 

- fig. 3 shows a block diagram summarising the operational phases of the 
discrimination procedure object of the present invention; 

30 - fig. 4 shows the diagram of a function enabling to obtain a fictitious parameter 
useful to the calculation of weights entering the spatial filtering made according to 
the procedure of the invention; 

- fig. 5 shows a curve of the merit parameter G n (6) (Array gain) versus the arrival 
direction on the azimuth plane, for a receiver employing a discrimination algorithm 

35 according to the known art; 
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- fig. 5' shows a table comparing the main merit parameters of the receiver scope of 
the present invention, in which it is included G n (0), with the corresponding 
parameters of receivers employing different discrimination algorithms according to 
the known art; 

- figures 6 and 6' represent a curve and a table, respectively relative to the same 
merit parameters of figures 5 and 5' calculated, contrarily to the first ones, 
employing reduced rank matrices; 

- fig. 7 shows a curve of the merit parameter G n (Q) versus the direction on the 
azimuth plane for the receiver employing the discrimination procedure object of the 
present invention; 

- fig. 7' shows a table making the comparison of the main merit parameters of the 
receiver scope of the present invention, with the corresponding parameters of 
receivers employing the discrimination algorithms according to the known art, said 
parameters being calculated using matrices differently reduced compared to those 
of figures 6 and 6'; and 

- fig. 8 shows a general block diagram of the receiver scope of the present invention. 
Detailed description 

Making reference to fig. 1, we notice an array antenna consisting of N 
electromagnetic field sensors identified a-i, a 2) a jt a N , of known type, arranged 
in straight line and separated one from the other by a distance d, typically d = /J/2, 
where k is the wave length of the radiofrequency in the middle of the band used by 
the particular mobile radio system employing the procedure scope of the invention (16 
cm approx. in case of GSM 900 MHz, 8 cm in case of DCS 1800 MHz). The array is 
invested by a given number of plane waves s 1t .... s„ s R% whose arrival directions 
on the azimuth plane 0i, ...,9/ are indicated for two of them, Values e n , 0/ 
corresponding to angles formed by radiuses s 1f ... t s, with the line of sensors a 1t a N . 
In a real scenario where the array is used by a base transceiver station (BTS) of a 
cellular telephone system, for instance GSM or DCS, plane waves s 1f .... s, 
correspond respectively to echoes of a wanted signal (identified by index / = u) 
transmitted on an assigned carrier and time slot, and to interfering echoes generated 
in neighbouring cells by communication channels which re-employ the same 
frequency in the same time slot. 

Fig. 2 is useful to evaluate the progressive delay, or the corresponding phase 
shifting, according to which a plane wave s,- invests the different sensors of the array. 
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If s(t) is a signal impinging on a uniform linear array from direction 9, the signal 
received by ihe n' h sensors x n (t) is a phase shifted replica of the s(t) according to the 
relation: 

,2{n-i)7idcosie) 

(D 



K with n = 1...N 

Thus the replicas on different sensors can be described by 



'(0 = 



1 





~<(e)~ 


5(0 = 


<(e) 




«C(e) 



s(<) = d-(eM/) 



where d* (9) is the response of the n th sensor to a signal coming from the direction 0. 

Defining the response d*(Q) as the complex conjugate of ^(9) is useful, in the 
following, to write linear combinations as scalar products. 

In case of more signals Sj(t) impinging from directions 0 M the signals x n (t) received by 
different sensors are described by: 



*„(/) 



= Z^(e l > / (/)- 



(2) 



where n n (t) are the uncorrelated noise components on array sensors 

Making reference to fig. 3, we notice four blocks in cascade called ARY, 

STDOA, STPESI, and FSPAZ, respectively. The ARY block is the sensor array of fig. 

1 and 2 which gives the signals x n (t) given by the expression (2). 

The STDOA block processes the signals x n (t) to estimate the number R of plane 

waves Si that invest the array ARY and the relative arrival directions (DOA Directions 

Of Arrival) on the azimuth plane, given by angles 9j. The operation of the STDOA 

block is known, for instance from [3] and [4] . 

The information obtained by the STDOA block, in particular the directions e it 

where implicitly a particular value / = u of index / indicates the direction 0 U of the 

wanted signal, is transferred to the next STPESI block that employs it in the 

calculation of appropriate spatial filtering coefficients w n , according to a novel method 

which shall be shortly described. 

Weights w n , together with the signals x n (t) transduced by the array ARY, are allowed 
to reach the FSPAZ block that gives a reception signal r(t) cleaned from interferents, 
thanks to a spatial filtering (beam forming) employing the following expression: 
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where it can be noticed thai the spatial filtering , is a simple linear combination of 
signals x n (t), transduced by the array ARY, weighed by coefficients w n . 

The (3) can be written in a form suitable to highlight the s, signals impinging on 
the sensors of the array ARY: 

( \ 

/ ^ n J n 

whose usefulness is that to enable the calculation of weights w n according to the 
procedure scope of the present invention. It must be pointed out that incident signals 
s,(t) become accessible only through the corresponding transduced ones x n (t).. 

For the determination of weights w n it is convenient to rewrite the (4) in vectorial 
form highlighting the contributions of the wanted signal s u (t), of interferents s*(t), and 
of the noise n(t), as it results from the following expression: 

/■(/) = s,(/)(d"w) + 2>,(/XcM + wn(/) (5) 

- where: simbol H indicates a transposed and conjugated matrix 

- w is the N element vector of coefficients w n ; 

- d H is the N element vector containing the response of the array in the direction 0 U 
of the wanted signal s u {t) ; 

- cf is the N element vector containing the response of the array in the i th interferent 
direction ; the index /' of the summation extends to value R-1;n(/) is the N 

elements noise vector. 

As for the operation of the STPESI block, it is convenient, with the aid of the (6), 
to introduce two gain parameters, respectively identified G n and G* in order to be able 
to introduce constraint conditions on said parameters for the calculation of weights w n 
. The above mentioned gains are the following: 
Gain of wanted signal on noise: 

r |d " w|r (6> 

is the gain in the S/N ratio versus the traditional case of use of a single element 
antenna. The notation || ... || indicates the Euclidean norm. The maximum value of GJs 
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equal to ||d| 2 

Gain of the wanted signal on the interferent /: 



G. = 



is the gain in the C/l ratio, concerning 



|d"w 



cfw 



(7) 



an interferent /', versus the traditional case of 
using a single element antenna, Analogously the gain of interferent on wanted signal 
is given by: 



J_ 
G. 



IcfwJ 



With these assumptions, the STPESI block calculates the vector of weights w 
satisfying the two following conditions: 

A) constrains the gain G„ in order that it is higher than, or equal to, a fraction a 2 of the 
maximum value, in order to obtain the above mentioned advantages. The 
parameter a 2 can be varied in a continuous fashion in order to trade off 
interference reduction for directivity towards the user The condition A) is equivalent 
to the expression: 



W 



■>* 2 lldf 



(8) 



w 



B) minimizes the sum of gains of interferents on wanted signal, maintaining the 
constraint A) to the purpose of reducing the interference in uplink and downlink. 
Condition B) is equivalent to the following expression: 



mm 



f 1 ^ 



= min 



(9) 



where C" is a matrix formed by vectors cf of the (5). 

Should all the interferents have the same power the (9) is equivalent to minimize the 
ratio between the sum of interferents" powers and the power of the wanted signal 
(//C). 

The. STPESI block calculates the vector of weights w imposing the conditions 
(8) and (9). To simplify the computation and without involving any conceptual limitation 
to the procedure, it is imposed that: 
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d"w = l (10) 
The above mentioned normalisation condition involves the multiplication of vector w 
by a constant scale factor, possibly complex, not influencing the spatial filtering made 
by the equation (5) since, as it can be noticed from the same, the wanted signal, 
interferents and noise would be multiplied by the same constant, and therefore the 
filtered signal, which can be denormaiized dividing it by the same scale factor. With 
the setting of the (10) the (8) becomes: 

1 



llwll 

and the (9) becomes: 



^>a 2 \\d\\ 2 (8') 



in(||c"wf) (9') 



10 mi 



The STPESI block employs the method of Lagrange multipliers to make a 
constraind minimization of (9'). To this purpose, it is first constructed a function 
F(wA,p,x) of variables w,/.,p,x . indicated with (1 1 ): 

( 1 ^ 

15 F(wA,P,x) = lC /f w| : + X flwjj 2 - — -r + 2p(Re(d H w)-l) + 2 X (lm(d /f w)-l) 

V a-\\d\\-J 

where parameters w,X,p,x are arbitrary constants to be determined equalling the 
gradient of F(w,X,p,x) to zero. The resulting equations are given by: 



— V K) = 2CC H w + 2Xw + 2d(p + Jx) = 0 



I dp> ' J d% J " K 1 

Solving for |3,x and replacing it is possible to obtain the value of vector w according to 
the X parameter, the resulting expression is the following: 

(cc'+xiV'd 

w = -^ (13) 

d"(CC" +U) d 

25 where I is a unit matrix of the same order of the matrix C. 



■WO 99/33141 




PCT/EP98/08165 



3 



To solve the (13) it is first necessary to aetermine the value of the parameter X. 
This real scalar is obtained assigning the desired value of the gain G„ given by the 
(8"), that is equivalent to remove the > sign from the same expression, having 
imposed a precise value of a 2 . The (8') therefore becomes: 

G " = jjf = or2 N 2 (8") 

Using a known factorization of eigenvalues and eigenvectors of the CC H 
matrix of the type: 

CC" = E*AE (14) 

where A is the diagonal matrix of the eigenvalues of CC", and E is the matrix of 
eigenvectors, we obtain the following expression of the gain (8"): 



~ = ™ (15) 



^(A,+xy 

where A, are the eigenvalues, real and not negative, of matrix CC" ; e, are the 
elements of the vector e = |e " dj more particularly, the elements e, correspond to the 
module of the components of vector E"d ; and X is the above mentioned real scalar, 
considered continuously variable to suit the arbitrary character of the selection of a 2 . 
From the (8") and (15) the following equation is reached in the unknown X : 



e: 



w 



whose solution gives the numeric value of parameter X that replaced in the (14) gives 
in its turn the vector of weights w. 

In view of a reduction of the computational complexity requested to the STPESI 
block, once the numeric value of the unknown X is obtained as said above, these can 
be replaced in the following equation: 



' 'WO 99/33141 




PCT/EP98/0816S 



_E(A + Uy i E f! d 
" "d w E(A + Xl)-'E H d (17) 
obtained replacing the (14) in the (13). The (17), contrarily to the (13), has the 
advantage not to require difficult matrix inversions, except for the easy inversion of the 
term in brackets which is a diagonal matrix. 

As far as the solution of equation (16) in the unknown value X. fig. 4 shows with 



a continuous line the diagram of the function (or equivalent^ cr (A) ) and 

Mi idii 



with a dotted line the diagram of the function to be maximised ! (\) that 

corresponds to the gain in C/l in case of equal power interferents. 

In fig. 4 the x-axis indicates the X value normalized at the value of the maximum 
) eigenvalue max(Aj) while on the y-axis values are in dB . 

As it can be noticed from fig.4 the function — is always increasing for X>0 

and tends to 1 for 1-kc 

The study of both functions suggests a procedure for the solution of the (16) in 
the unknown X . To this purpose the following steps are foreseen: 

- we calculate — i--(x = o); 

H Ml 

' f |w| 2 iidj| 2 ^ = °^ " a2th6n WS choose x= °- because it is already the desired value 
since for X = 0 we obtain the unconstrained minimum value of |c"w| 2 ; 

' f Iwjpildfl 2 ^ = °^ * a2then X is incre ased from 0 employing an iterative method for 

equation solution. In particular using the first order Newton method we obtain the 
iterative relation: 
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f 

df 



f 1 



where 



function f(x) is properly selected in order to guarantee the sequence 



x- 1 



convergence. Practically f(x) can be chosen as f(x) = log(x) or f(x) = 

x + 0.05 

For example from the diagrams of fig. 4 it results that the case represented 
corresponds to a very low noise gain G„ for X = 0. On the contrary, selecting a 1 = 0.5 
(-3dB) so that the gain in S/N is only 3 dB below the maximum achievable value, it 
results that the solution of (16) is approximately X = 0,15xmax(A). 

Referring to figures 5, 5\ 6, 6', 7 and T the performance of the discrimination 
procedure scope of the present invention are now compared with similar 
performances of the main beamforming algorithms known so as described in [2]. The 
scope is that to highlight additional drawbacks of the known art now superseded by 
the present invention. 

Merit parameters, expressed in dB, highlighted in the tables of figures 5', 6' and 
T attached to the relative figures 5, 6 and 7 are the following: 

- the gain G n in the S/N ratio, versus a traditional antenna having a single element, 
given by (6); 

C 

- the gain in the - ratio between the wanted signal and the interfered signal at the 

output of the beamformer (FSPAZ block), versus a traditional antenna having a 
single element; 

.... S 

- the gam in the - - ratio between the wanted signal and the noise plus 

the interfered signal at the output of the beamformer (FSPAZ block), versus a 
traditional antenna having a single element. Concerning this last parameter, the 
following expression applies: 

S oui = S,„ 
+/ M. — /,._ (21) 
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Using the same terminology of [2] , the known algorithms considered are the 

fnllnwino* 

1. Minimum Variance (MV): known the wanted signal direction, the weights w n are 
selected in order to minimize the power at the beamformer output with the 

5 condition d w = 1 . 

2. Linear Constrained Minimum Variance (LCMV): in addition to the condition 

d H w = 1 , the weights w n minimize the power at output of the beamformer with the 
additional constraints that the array response is null in the interferents directions 
c"w = 0. 

10 3. Linear Constrained (LC): weights w n meet the previous constraints on the wanted 
signal and on interferents, also minimizing ||w|| : . 

The procedure scope of the present invention is indicated in the tables with: 
4. Constrained Gain Minimum Interference (CGMI): weights w n minimize the sum of 

gains of interferents on wanted signal (jjc" wf ) f with the constraint that the gain G n 

15 of the wanted signal on the noise is higher than a pre-set threshold (6 dB in the 
following examples). 

To obtain the above mentioned tables and comparison diagrams the following 
scenario was selected, particularly suitable to highlight the defects of the different 
algorithms: 
20 • an array of 8 elements is used 

• a wanted signal with arrival direction orthogonal to the array plane (DOA = 0°); 

• two interferents of equal power, one half of the wanted signal , with a small angular 
separation from this last: (C/l = 0 dB; DOA = [-5°, 5°]); 

• very low noise level (S/N = 40 dB); 

25 • arrival directions of the wanted signal and interferents perfectly known. 

Making reference to the comparison table 1 of fig. 5\ we can notice that the first 
3 algorithms show poor performance due to a really low gain G n (-25 dB for LCMV and 
MV and -5 dB for LC), contrarily to what takes place for the CGMI method, which has 
a constrained G n equal to +6 dB ( used in the discrimination procedure scope of the 

30 present invention. This means that, when using the LCMV or MV algorithm, the power 
radiated in downlink towards the mobile is 25 dB lower compared to the case of single 
element antenna. In this condition, also if interference is absent, the S/N at the mobile 
could be too low for a correct detection. 
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In the corresponding fig. 5, showing the radiation diagram that is the gain c (9) 
cf the array concerning the application of the algorithm LCMV, we notice that the 
above mentioned diagram shows a very low value (-25 dB) in the wanted signal 
direction, null (-co dB) in the directions of interferents, and high values in the other 
5 values of 0. Thus fig. 5 shows how, with algorithms from known art, most of the power 
is radiated in directions different from the wanted direction. 

This behaviour is due to the fact that the LCMV algorithm, as the other known 
algorithms can run up against the inversion of almost singular matrixes, that is with 
very high condition number, when arrival directions of interferent and wanted signal 
io are close each other. 

A way to mitigate this effect is to use matrixes of reduced rank, obtained 
selecting only the highest singular values of the matrixes to invert. 

Figures 6 and 6', maintaining the previous scenario, show the effects obtained 
employing reduced matrixes in the different algorithms indicated; the reduction is 
15 conducted selecting only the 2 highest singular values. Concerning the CGMI method 
used in the invention, it does not involve the inversion of nearly singular matrixes, 

The radiation diagram of fig. 6, still showing the array radiation diagram g„(6) 
concerning the application of the algorithm LCMV, has a much more higher directivity 
in the wanted, however, also interferent signals result very high. Table 2 of fig.6' 
20 confirms this behaviour and shows that with the known algorithms a high g„ value is 
obtained in the wanted signal direction, but the capability to null the interferents is lost, 
as it results from the low values of parameter A(C/I). On the contrary, with the CGMI 
method, in which the gain c m is previously set at 6 dB, a higher reduction of 
interferents is obtained. 
25 Figures 7 and 7\ maintaining the previous scenario, show the effects obtained 

increasing to 3 the number of the highest singular values used in the inversion, in 
order to partially recover the capacity to nullify the interferents, though accepting to 
lose something in the gain g„ since we come mgre close to the situation of figures 5 
and 5'. Actually, the results of Table 3 of fig. T show the correctness of the reasoning, 
30 but show also that the gain G n of known algorithms is too low for a correct antenna 
operation, contrarily to what takes place for the CGMI method used in the invention. 

Fig. 7 shows the radiation diagram obtained applying the CGMI method used in 
the present invention, confirming the results shown in the table. 
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What we can conclude from all comparisons made is that, differently form 
known algorithms, the algorithm which is the scope of the present invention allows to 
easily and continuously trade off antenna gain in the wanted signal direction for 
capacity to null the interferents. The conclusion is supported by the fact that the 
5 passage from the situation of Table 2 (fig. 6') to that of Table 3 (fig. 7'), describing 
known art results, does not allow continuously variable intermediate situations, since 
reduced matrices can be only obtained selecting an integer number of singular values. 

Fig. 8 gives a general representation of a receiver for base transceiver station of 
a cellular telecommunication system GSM, or DCS type, employing the discrimination 
10 procedure scope of the present invention. More particularly, the receiver of fig.8, and 
the relative associated transmitter, can find application in frequency division (FDMA) 
or time division (TDMA), or mixed FDMArTDMA multiple access systems re-employing 
a same frequency group in adjacent territorial areas. 

Making reference to fig.8, we can notice the array antenna ARY consisting of N 
15 elements a-,, .... a jf .... a Nl connected to a block RIC-FDMA/TDMA from which digital 
demodulated signals CHrrij come out at each time slot, reaching a process module 
PROC. In the symbol CHrrij, the index m indicates the generic carrier m-th of the M 
separate carriers assigned to the receiver, while the index j indicates the j-th of the N 
replica of the reception signal relative to the m-th carrier. A snapshot is thus supplied 
20 of the signals coming out from the above mentioned block at the current time slot. The 
process module PROC processes, according to the procedure described above, the N 
replica of each one of the M channels CHm and supplies M digital signals CHFm, 
spatially filtered, corresponding to the M channels simultaneously received. 

In the operation, the RIC-FDMA/TDMA block performs all the operations 
25 necessary to the reception of M channels of the FDMA/TDMA type from each one of 
the N sensors of the array, that is: 

• radiofrequency filtering for the suppression of the spurious out of the total reception 
band and subsequent frequency splitting; 

• conversion at intermediate frequency and filtering; 
30 • conversion analogue to digital; 

• best demodulation and filtering to obtain demultiplexed signals FDMA/TDMA, to be 
sent to the PROC block. 

The M channels CHF1 , CHFm, .... CHFM, spatially filtered, coming out from 
this last block suffer the following additional processing inside the receiver of fig.8: 
35 • reconstruction in base band of the M original transmission bursts associated to the 
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M channels and positioning inside frames and multiframes; and 
• new coding of digital bursts in a format (PCM 30 channels) compatible with the 
protocol (LAPD) adopted on beams connecting the base transceiver stations to the 
relative station controller (BSC). 
5 The original character of this receiver is of course that of the PROC block, so 

that the RIC-FDMA/TDMA block can be considered known to the field technician. 
More particularly, it is possible to adopt for the above mentioned block, an architecture 
based on single narrow band receivers, as for instance results being that of the 
disposition of antenna according to the known art mentioned above, or a disposition 
10 based on the use of broad band transceivers simultaneously processing several 
channels, like for instance the one described in the European patent application No. 
97830229.7 filed under the name of the same applicant. 

Profitably, the PROC block of fig. 8 can be implemented through a 
microprocessor for the mathematical processing of digital signals (DSP), or more 
15 adequately, through digital integrated circuits of the ASIC type (Application Specific 
Integrated Circuit). 

Therefore, while a particular embodiment of the present invention has been shown 
and described, it should be understood that the present invention is not limited thereto 
since other embodiments may be made by those skilled in the art without departing 
20 from the scope thereof. It is thus contemplated that the present invention 
encompasses any and all such embodiments covered by the following claims. 
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CLAIMS 

1. Discrimination procedure of wanted signal (5- u (/)) from a plurality of 
cochannel interferents (s t (i)) received by array antennas (ARY) of frequency division 
or time division or mixed multiple access telecommunication systems, re-employing a 

5 same group of frequencies in adjacent territorial areas, including an estimate phase 
(STDOA) of the number of said interferents and of the their arrival directions 9 f and of 
the arrival direction of said wanted signal 0 Ut and a subsequent spatial filtering phase 
(FSPAZ, STPESI) in which signals x n (t) transduced by relative sensors (a 1( .... a N ) of a 
said array (ARY) are linearly combined among them through multiplication coefficients 
10 w n , or weights, obtaining a reception signal r(t) cleaned from interferents; said phases 
being repeated for each one of the time slots, in the TDMA frame, in which the wanted 
user transmits, characterized in that said weights meet the following conditions: 

A) the gain of said array in the ratio between said wanted signal and the noise 
components after spatial filtering, compared to the use of a said single sensor, is 

15 constrained in order to result higher than or equal to a properly selected threshold; 

B) the sum of gains of each said interferent on the wanted signal after spatial filtering, 
compared to the use of a said single sensor, is minimized . 

2. Procedure according to claim 1, characterized in that said spatial filtering is 
made employing the following expression valid for complex envelopes in base band: 

20 r(f) = ^„x n (t) 

n 

where: 

- r(t) is a filtered signal whose characteristics reflect said discrimination of the wanted 
signal from said interferents; 

- x n (t) are the said signals transduced by the relative said sensors {a, a N )\ and 

25 - w n are said multiplication coefficients, or weights. 

3. Procedure according to claim 2, characterized in that said expression used to 
perform the spatial filtering is equivalent to the following one, given in vectorial form: 

r{t) = s u (t)(d H w) + X s s (/Xcf w) + wn(/) 

where: 

30 - s u (t) is said wanted signal inciding on said array (ARY); 

- Si(t) is an interferent Mh of said plurality inciding on said array (ARY); 

- w is the N elements vector of said weights w n , relative to said sensors (ai a N )\ 
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- d H is the /V-element vector that indicates the response of the array in the direction 
of the wanted signal; 

- c" is athe A/-eiement vector that indicates the response of the array in the 
interferent direction /-th of said plurality; and 

- n(/) is a noise vector at N elements, each one associated to one said sensor. 

4. Procedure according to claim 3, characterized in that said condition A) is 
equivalent to the expression: 



lld"w 



>a 2 \\d\\ 



l|W|| 

and said condition B) is equivalent to the expression: 



'"C'w"^ 



mm 



Id"wn j 



where: 



- a' is said pre-set fraction of said maximum value |d|| 2 = max 



d WI; 



being 



. ||...|| the Euclidean norm; and 

- C H is a matrix formed by vectors cf . 

5. Procedure according to claim 4, characterized in that said spatial filtering is 
made imposing the following normalization constraint: 

d*w = l. 

6. Procedure according to claim 5, characterized in that said vector of weights w 
is calculated through the following expression: 

(CC* + uY& 

w "d /, (cc H +xi)" 1 d 

where: 

- I is a unit matrix of the same order of said matrix C ; 

- X is a real scalar parameter obtained solving the expression related to said 
condition A), in presence of said normalization constraint and imposing a precise 
value to said pre-set fraction a 2 . 

7. Procedure according to claim 6, characterized in that said matrix CC H is 
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- A is the diagonal matrix containing the eigenvalues of said matrix , and 

- E is the matrix of the eigenvectors. 

8. Procedure according to claim 7, characterized in that said vector of weights w 
is also calculated through the following expression: 

= E (A + M)~'E*d 
W " d w E(A + JJ)"'E w d ' 

9. Procedure according to claim 7, or 8, characterized in that said parameter X 
is obtained solving the following equation having one unknown value: 



- A, are the eigenvalues, real and not negative, of the matrix CC" ; and 



10. Procedure according to claim 9, characterized in that the above mentioned 
equation in said unknown X is solved in iterative way through the following steps : 




where: 



- e i are the elements of a vector e = E w dl 



- it is calculated 



iwfifdir 



;tU = o); 




(A = 0) > a 2 then we choose A. = 0, because it is already the desired 



value since from the per X = 0 we obtain the minimum value of ||c 




20 



- if 



■j (a = 0)<a 2 , X is increased from 0 using the iterative relation: 
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0 



df 



f 1 



: a — w 

where function f(x) is properly selected in order to guarantee the sequence 
convergence 

Iterations are executed until it results — I a. 2 

iwfiidi 2 " 

11. Procedure according to any of the previous claims, characterized in that 
N = 8 is said number of sensors (a,. a«), and a 2 = 0.5 is said pre-set fraction of 
said maximum value ||d|| 2 . 

12. Procedure according to any of the previous claims, characterized in that 
said telecommunication system is a cellular telephone system. 

13. Receiver for frequency division, or time division, or mixed multiple access 
telecommunication systems, which re-employ a same frequency group in adjacent 
territorial areas, applying the discrimination procedure object of the previous claims, 
including: 

- an array of sensors of electromagnetic field (ARY); 

- radiofrequency filtering means (RIC-FDMA/TDMA) of the signal coming from each 
sensor (a 1f a w ) for the suppression of the spurious out of the global reception 
band, having downstream splitting means of a plurality of carriers; 

- means for conversion and filtering (RIC-FDMA/TDMA) at intermediate frequency; 
analogue-to-digital conversion means; demodulation and filtering means to obtain 
in base band demultiplexed signals FDMA/TDMA of relative communication 
channels; 

- process means (PROC) of said signals in base band and of reconstruction of the 
information originally transmitted on said channels; 

characterized in that said processing means include also: 

- estimate means (STDOA) of the arrival direction 6 U of a wanted signal s„(t), and of 
the number and arrival directions 9j of a plurality of cochannel interferents signals 
Si(t); 

- spatial filtering means (FSPAZ) of said signals in base band, to obtain a reception 
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signal r(t) cleaned from said interferents, said means linearly combining among 
them x n (t) signals corning from the relative said sensors (a 1t a N ) of a said array 
(ARY), through multiplication coefficients w n , or weights; 

- calculation means (STPESI) of the gain of the wanted signal on the noise, spatially 
5 filtered, compared to the traditional case of utilization of a single element antenna; 

- calculation means (STPESI) of the gain of the wanted signal on the interferent /-th, 
spatially filtered, relating to an interferent /-th and for all the interferents, compared 
to a traditional case of use of an antenna to a single element; 

- means calculating said weights (STPESI) imposing a value of said gain of the 
10 wanted signal on noise, higher than or equal to a duly selected threshold, and 

simultaneously a minimum value of the sum of gains of each said interferent on the 
wanted signal. 

14. Receiver according to claim 13, characterized in that N = 8 is said number of 
sensors (ai, a*), and a 2 = 0.5 is said predetermined fraction of said maximum 

15 value. 

15. Receiver according to claim 13, or 14, characterized in that it is used in a 
cellular telecommunication system. 

16. Receiver according to any claim 13 to 15, characterized in that said process 
means (PROC) consist of digital integrated circuits of the ASIC type {Application 

20 Specific Integrated Circuit). 




Fig. 2 



WO 99/33141 PCT/EP98/08165 

2/6 



ARY 




STDOA 


6i 


STPESI 


W; 

J ™> 


FSPAZ 






-> 





r(t) 



Fig. 3 



■ 


1 1 


n 1 1 

vc/i ; 






! /Va 2 










.7 (dB) • 
















.1 


I 1 '• 

[ " 


—I — I i 1 







-0.5 



0 0.15 0.5 1 
X /max(A) 



1.5 



Fig. 4 



SUBSTITUTE SHEET (RULE 26) 



WO 99/33141 




PCT/EP98/08165 



io r G n(Q) 




e 



Fig. 5 





Gn 


A(S/N) 


A(Cl\) 


A(S/(I+N)) 


LCMV 


-25.3 


-24.8 


232.9 


15.1 


LC 


-5.3 


-5.8 


289.8 


34.1 


MV 


-25.9 


-25.5 


17.9 


13.4 


CGMI 


6.0 


5.2 


3.5 


3.5 



Fig. 5' 



■ " WO 99/33141 




PCT/EP98/08165 




Fig. 6' 

SUBSTITUTE SHEET (RULE 26) 




-80 -60 -40 -20 0 20 40 60 80 

deg ree 



Fig. 7 





Gn 


A(S/N) 


A(C/I) 


A(S/(I+N)) 


LCMV 


-5.3 


-5.0 


296.3 


35.1 


LC 


-5.3 


-5.0 


289.9 


35.1 


MV 


-6.1 


-5.8 


14.6 


14.6 


CGMI 


6.0 


6.5 


3.7 


3.7 



Fig. T 



WO 99/33141 PCT/EP98/08165 

6/6 



LL 
O 



E 

Li. 
O 



LL. 
X 

o 




A A A 



o 



o 



o 















e" 


E 


E 


X 


X 


X 


o 


o 


o 



o 



















o 


o 



00 

d 

MB 



or 



CO 



CM 
CO 



> 



> 



co 



co 



> 



INTERNA 



NAL SEARCH REPORT 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 6 H01Q3/26 



According to International Patent Classification i\PC) or to both national Hac 



Int .tional Application No 

PCT/EP 98/08165 



B. FIELDS SEARCHED 



Minimum documentation searched {classification system followed by classification symbols! 

IPC 6 H01Q 



Documentation searched other than minimum Documentation to the extent that such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and, wnere practical, search terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category J Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



EP 0 809.323 A (CHOI) 26 November 1997 

see page 12, line 44 - page 13, line 23; 
figure 7 

see claims 1-31; figures 1-6B 

WO 95 22873 A (ERICSSON) 24 August 1995 

see page 3, line 14 - page 4, line 25; 
claims 1-52; figures 1-5B 

FR 2 747 792 A (LE CENTRE THOMSON D APPLIC 

RAD) 24 October 1997 

see abstract; claims 1-7; figures 1-3 

US 5 574 824 A (SLYH ET AL. ) 

12 November 1996 

see abstract; figures 1-5 

-/-- 



I- 3, 

II- 16 



I- 3. 

II- 16 



1.16 



1.16 



m 



Further documents are listed in the continuation of box C. 



Patent family memoers are listed in annex. 



° Special categories of cited documents : 

"A" document defining the general state of the art wnich is not 
considered to be of particular relevance 

"E" earlier document but published on or after the international 
filing date 

"L" document which may throw doubts on priority claim(s) or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

"O" document referring to an oral disclosure, use. exhibition or 
other means 

"P" document published prior to the international filing date but 
later than the priority date claimed 



"T" later document published after the international filing date 
or priority date and not in conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

"X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y" document of particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other sucn docu- 
ments, such combination being obvious to a person skilled 
m the art. 

"&" document member of the same patent family 



Date of the actual completion of the international search 

21 April 1999 



Name and mailing address of the ISA 

European Patent Office. P.B. 5818 Patentlaan2 
NL - 2280 HV Rijswijk 
Tel. (+31-70) 340-2040. Tx. 31 651 epo nl. 
Fax: (+31-70) 340-3016 



Form PCT/ISA/210 (second sheet) (July 1992) 



Oate of mailing of the international search report 



28/04/1999 



Authorized officer 



Angrabeit, F 



oaqe 1 of 2 



INTE 



IONAL SEARCH REPORT 



C.(Continuation) OOCUMENTS CONSIDERED TO BE RELEVANT 



Inl .tional Apptication No 

PCT/EP 98/08165 



Category * Citation ot document, with indication, where appropnate. of the relevant passages 



EP 0 604 956 A (NEC) 6 July 1994 
see abstract; figures 1— I I 



Relevant to claim No. 



Form PCT/ISA/210 (continuation of second sheet) (Jufy 1992) 



page 2 of 



2 



INTERJ^VIONAL SEARCH REPORT 

Information on patent family members 



tntt^^;tonal Application No 

PCT/EP 98/08165 



Patent document 
cited m search report 



Publication 
date 



Patent family 
rnember(s) 



Publication 
date 



EP 


0809323 


A 


26- 


-11-1997 


CN 


1171664 


A 


28-qi -i qqq 












JP 


10098325 

-I \J \j Z* {J *J C— ~J 


A 


— • V ~T i J JO 












US 


5808913 


A 


15-OQ-IQQft 

J. W 1 17 JO 


WO 


9522873 


A 


24- 


-08-1995 


US 


5566209 


A 


15-10-1996 












AU 


694175 


B 


16-07-1998 












AU 


1723795 


A 


. 04-09-1995 












CA 


2159712 


A 


24-08-1995 












CN 


1124548 


A 


12-06-1996 












FP 




A 


9 A— ft 1 — 1 OOR 












FT 




A 


1U Ijjj 












JP 


851 1 ^Q4 


T 

i 














SG 


48507 


A 


17-04-1998 


FR 


2747792 


A 


24- 


•10-1997 


NONE 








US 


5574824 


A 


12- 


■11-1996 


NONE 








EP 


0604956 


A 


06- 


07-1994 


JP 


2663820 


B 


15-10-1997 












JP 


6204902 


A 


22-07-1994 












US 


5689528 


A 


18-11-1997 



Form PCT/ISA/210 (patent family annex) (Juiy 1992) 



